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1 Introduction 



A full understanding of the properties of a quantum system and prediction of its behaviour 
cannot be achieved without knowing its spectrum, i.e. the energies of its bound states 
and resonances. The problems of this kind emerge not only in fundamental research 
concerning particles, nuclei, and atoms, but also in engineering. For example, modern 
semiconductor devices based on nano-structures, cannot be properly designed without 
accurate treatment of various transitions among stationary and quasi-stationary states of 
the charge carriers (see, for example, Ref. [T| and references therein). 

The problem of locating bound states is as old as quantum mechanics itself. Over 
a century, plenty of exact and approximate methods were developed for solving it. The 
notion of quasi-stationary states (or quantum resonances) emerged at later stages of the 
development of quantum mechanics. The attention to such states was drawn by Gamow 
in his pioneering works on the a-decay [2]. 

At the beginning, the development of the methods for locating resonances was hin- 
dered by computational difficulties, which are more challenging than that of the bound 
state problem. The progress in this field was therefore delayed till the advent of modern 
computers. Another consequence of these difficulties was that the methods for solving the 
bound and quasi-bound state problems were developed separately despite the fact that 
the bound and resonant states have essentially the same mathematical nature. Indeed, 
they correspond to the S-matrix poles on the complex energy surface at all of which the 
asymptotics of the solutions of the Schrodinger equation has the same functional form, 
namely, the outgoing wave. The only difference between bound and quasi-bound states 
is that the corresponding poles are situated in different domains of the energy surface. 
As a result, the wave function asymptotics look deceptively different. Being found at 
real negative energies, the bound state solutions at large distances exponentially decay. 
However, the decaying exponential function is the same outgoing wave, but taken with 
pure imaginary momentum that corresponds to a negative energy. 

Another fact that makes the kinship between the bound and quasi-bound states appar- 
ent, is their mutual transformations. This happens when one increases or decreases depth 
of the potential. A gradually deepening potential "sucks" the resonances in. The corre- 
sponding poles of the S'-matrix move towards the threshold point {E = 0) and eventually 
cross it over to the bound state domain. 

Common mathematical ground of the bound and resonant states implies that it should 
be a unified way of locating them. In quest for such a way most of the modern methods 
were developed. A review of the existing methods for solving the quasi-bound state 
problem can be found, for example, in Ref. [3]. One of the new unified approaches not 
mentioned in that book, is based on a direct calculation of the Jost-function at complex 
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energies [1] (for references concerning recent development of this method, see also Ref.[S]). 

All the methods for locating resonances can be divided in two groups. One of them 
comprises the approaches, in which all the calculations are done only at real energies. 
Such methods are rather simple but have limited abilities. They usually fail for wide as 
well as extremely narrow resonances and are not unified in the abovementioned sense. The 
other group combines the methods based on locating the S-matrix poles (or Jost function 
zeros) on the complex energy surface. These approaches are very powerful, accurate, 
and unified. The price one has to pay for that is that they require rather sophisticated 
calculations. 

Even the unified complex-energy methods are not fully universal. They are not appli- 
cable, for example, to potentials that cannot be analytically continued to complex values 
of the distance r. They have difficulties in dealing with energy- dependent and non-local 
potentials. The methods based on the complex rotation, cannot reach the so-called vir- 
tual states that correspond to the S-matrix poles at negative imaginary momenta, i.e. at 
negative energies on the unphysical sheet of complex energy surface. 

Thus, one can choose either a simple but inaccurate real-energy method or a powerful 
but complicated complex-energy one. None of them, however, is universal. It would be 
desirable to have a method that could combine in itself both the simplicity of the real- 
energy approaches and the power and accuracy of the complex-energy ones. In this paper, 
we suggest such a method. Besides this, the proposed method is insensitive to the nature 
of the interaction potential and therefore can be considered as universal. 

The idea is based on the fact that two analytic functions coinciding on a curve segment, 
are identical everywhere in the complex plane (the so-called coincidence principle^). 
Therefore if we calculate (using any appropriate or available method) the S'-matrix along 
the real positive axis of the energy plane and accurately fit the obtained values with a 
meromorphic function, for example, with a Pade function, then we can expect that this 
approximate function will have practically the same singularities at complex energies as 
the exact S-matrix. In this way, we can locate the spectral points (bound, virtual, and 
resonance states) by locating these singularities. In the suggested method, all the poles of 
the Pade function are located at once, because their coordinates are the fitting parameters 
and thus they are determined within the fitting procedure. 

By constructing the Pade approximant, we actually do an analytic continuation of 
the S'-matrix from a segment of the real axis, where it is given at a number of discrete 
points, to complex energies. Similar (but different) procedures in which function values 
on the complex plane are obtained using the knowledge of its values at a discrete set of 
real points, were more than once used before and proved to be successful. As an example, 
we can mention Ref. |7j where a potential V{r) numerically calculated at a set of real 
points ri, r2, . . . , rjv was used within a complex rotation method for locating resonances. 
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Another example is the Pade approximation of the Titchmarsh-Weyl m-function that was 
used in Ref. [S] also for locating resonances. 



2 Rational approximation of the /S-matrix 

Let us assume that we know the complex- valued function Si{k) for all real values of the 
collision momentum k G [0,oo). Alternatively, we can have the real phase-shift function 
6i{k) related to the S'- matrix as 



In practical calculations the interval [0, oo) is reduced, of course, to a finite segment 
[^min, ^max] with Sufficiently small fcmin and large /cmax, that covers all significant oscillations 
of the function Si{k). 

We are going to find an approximate function Si{k) such that the difference Si{k) — 
Si{k) is minimal on the segment [/cmin, ^max]- Keeping in mind the well-known fact (see, 
for example, Ref. [6]) that two analytic functions are identical everywhere if they coincide 
on a continuous segment, we then expect that the approximate S'-matrix Se{k) will have 
almost the same singularities (resonance and bound state poles) as the exact one. The 
more accurately we approximate the S'-matrix on the real axis, the less different will be 
the poles of Si{k) from the corresponding poles of Si{k). 

There are many ways of approximating Se{k). A choice for the functional form of 
Se{k) is determined by the fact that it should have simple poles, i.e. at certain points fcj 
be proportional to ~ (A; — ki)~^ . Such behaviour is provided by a ratio of two polynomials 



which in numerical analysis is known as the Pade approximation of the order [N, M] when 
the parameters am and 6„ are chosen to reproduce the exact function and all its derivatives 
up to the order M + at A; = (see, for example, Refs. [HI [TO]). 

Here, we will use a special form of the approximation ([2]). Firstly, we observe that 
at high energies the S'-matrix tends to unity [llj, which can only be achieved if both 
the numerator and denominator polynomials have the same degree, i.e. M = N, and 
O'M = ^N- Secondly, at zero energy the S'-matrix is unity as well [TTj, which implies that 
Oo = &o (without loosing the generality, we can assume that Oq = &o = 1). Thirdly, the 
fitting parameters cannot in our case be found using the S'-matrix derivatives, which are 
not available. Instead, we use the algorithm described next. 



S,{k) = exp \i25,{k)] . 



(1) 



Oo + aik + a2k'^ + . . . + a^fk 



(2) 



6o + hik + h2k^ + ... + h^k^ 
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3 Fitting parameters 



Sdk) = , (4) 



To begin with, we re-arrange the polynomials in the form that includes their zeros explic- 
itly, namely, 

^ k-a 

SdE) = Ut-^- (3) 

n=l Pn 

Knowing some general properties of the S*- matrix, we can simplify the expression ([3]). This 
can be done as follows. First of all, we notice that the exact ^-matrix is the following 
ratio 

where f^^^^"^^\k) are the Jost functions |12j that are the amplitudes of the incoming and 
outgoing waves in the asymptotics or the radial wave function, 

Mk,r) ;^ft\k)h\~\kr) + ft'''\kW\kr) , (5) 

with the incoming and outgoing spherical waves being represented by the corresponding 
Riccati-Hankel functions h\~\kr) and h^^\kr). 

Now, we use the fact that fi'^\k) and f[°^^\k) are not independent. Indeed, the 
incoming and outgoing waves swap their roles when the momentum k changes its sign 
and also under the operation of complex conjugation. It is not difficult to show (see, for 
example, Ref. [12]) that this implies the symmetry properties summarized in Fig. [H In 
particular, we have 

fr'^{k) = ft\-k), (6) 

and thus 

This means that not all parameters in Eq. ([3]) are independent. The numerator polynomial 
should be the same as the denominator one but taken with negative k. Therefore 

~s.{k) = {-If n 1^ • (8) 

n=l Pn 

As compared to Eq. (j3]), this not only reduces the number of fitting parameters in half, 
but also improves the quality of the approximation since the correct structure of the 
5'-matrix, given by Eq. (JTj), is taken into account. 



5 



The procedure of finding the parameters /3„ consists of two stages. At the first stage, 
we once again re-arrange the polynomials in Eq. (|H]) to the form ([2]), 

N 

1 + J2 

Seik) = , (9) 

l + ^(-l)"a„r 

n=l 

where the numerator and denominator are divided by the product /3i/32 ■ ■ ■ Pn, which gives 
o-o = bo = I. In order to find the parameters a„, we multiply Eq. by the denominator 
of its right hand side and re-write it as 

N 

+ i-ir+'S,ik)] k^an = S,{k) - 1 . (10) 

n=l 

The last equation, taken at different points ki, k2, . . . , k^, constitutes a system of linear 
equations 

N 

'^Amnan = Bm, 771 = 1, 2,..., N (11) 

n=l 

determining the unknown coefficients a^. The matrices of this system are defined as 



A 



1 + (-l)"+i5,(U k:^ (12) 



and 

5„ = S,ik^) - 1 (13) 

with Se{km) being known values of the ^-matrix on the interval [fcmin, ^max]- Therefore, 
by solving matrix equation ffTTj) . we can determine the parameters a„. 

What we actually need are the S'-matrix poles, i.e. the parameters So, at the 
second stage, the parameters /3„ are determined as the complex roots of the polynomial 

N 

PN{k) = l + Y^i-ir ar^k^ . (14) 

n=l 

There are many robust and fast algorithms for finding such roots (see, for example, Ref. 
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4 Meaningful and spurious poles 



What are the mathematical and physical meaning of the poles of the approximate S- 
matrix Si{k)l As we know, all the poles of the exact S'-matrix are the spectral points, 
i.e. the bound and virtual staes (if any) as well as the resonances. Is this also valid for 
the poles of Se{k)7 

At a first sight, it seems that we should give a positive answer to this question. Indeed, 
the functions Se{k) and Sg{k) coincide on a segment of the real axis and thus, according 
to the coincidence principle of the complex analysis, they should be identical everywhere 
on the complex fc-plane. There are, however, two facts that cast some doubts on this 
reasoning. 

Firstly, the function Si{k) coincides with Si{k) not on a continuous segment, but only 
at a finite number of discrete points ki, k2, ■ ■ ■ , kjsf- And secondly, the number N of these 
points and thus the number of poles that the function ([8]) may have, is chosen arbitrarily. 

It is natural to expect that the smaller is the difference between Se{k) and Si{k) on the 
chosen segment, the better is the approximation of the S'-matrix at complex momenta. 
This difference can be made smaller by increasing the number of the points, in other 
words, the number of poles of Se{k). Surely, not all of these poles have physical meaning. 
Some of them may appear at "wrong" places, because of the non-zero difference S£{k) — 
Se{k) in between the points ki, k2, ■ ■ ■ , /cat. On the other hand, some of them must be close 
to the "true" poles, otherwise the approximate S-matrix would be too different from the 
exact one. 

Those poles of S^^k) that are close to the "true" S'-matrix poles are meaningful, while 
all the other poles are spurious. It is obvious that the positions of the meaningful poles 
are "tied" to the the "true" poles, while the spurious poles are "free to move" and appear 
randomly depending on the choice of the fitting points ki, k2, ■ ■ ■ , k^- 

Therefore, there is a simple way to distinguish meaningful poles from the spurious 
ones. Indeed, repeating the calculations with different number of fitting points, we can 
easily find those poles that appear more or less at the same positions. All the other poles 
should be regarded as spurious. In this way, we can also find an appropriate number A^, 
with which the meaningful poles converge to the "true" poles within a required accuracy. 

When doing the calculations, we found (see the next section) that the spurious poles 
in many cases appear in symmetric pairs that almost exactly cancel each other in the 
product (jH]). For example, if a spurious pole appears at = then it is accompanied by 
another pole at /c = /3m such that ~ —f^n- As a result, in the product (jH]), we have the 
f cLctor 

(fc-/3„)(fc-/3m) ^ ' ^ ^ 
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which is practically unity everywhere except the immediate vicinity of the spurious poles. 



5 Numerical examples 

The basic idea of the proposed method rests on a rigorous mathematical fact of the 
identity of two functions that coincide on a continuous curve segment. In numerical 
calculations, however, we can only guarantee that the exact and approximate functions 
Se{k) and S^lk) coincide at a number of discrete points along the real fc-axis. Of course, 
after the parameters of Se{k) are fixed, we can always check how this approximate function 
reproduces the exact one at the intermediate points. Our hope is based on the well-known 
fact that rational interpolation of the type works very well even with just few matching 
points|^. Certainly, the accuracy of the approximation is best near the matching points 
and should deteriorate when we move far away into the complex plane. 

In order to test how this rational interpolation works in our method, we performed 
calculations for several well-studied potentials, whose spectral points are known. Since 
in the proposed method nothing special is associated with the angular momentum, we 
tested here only the case i = 0. The exact values of the S'-matrix at the fitting points 
were calculated using the Jost function method described in Ref. The same method 
was used to locate the the exact spectral points, with which the approximate S'-matrix 
poles were compared. 

The first of the testing potentials is an exponential well. 



with ly(^) = 154.06 MeV and Ri = 0.76 fm, which roughly describes the S'-wave proton- 
neutron interaction in the triplet state (total spin = 1). With h'^ / (2m) = 41.47 MeV • fm^ 
this potential generates a bound state (the deuteron) at E = —2.2244674752 MeV. With 
£ = 0, it does not generate any other poles of the S'-matrix (at least within any physically 
reasonable domain of the complex plane). 

Since there is only one "true" pole, we expect that it is not necessary to have many 
fitting points in order to reproduce this pole using the Fade approximation. And indeed, 
starting with just one point, we see a very rapid convergence, so that at = 10 the exact 
value of the binding energy is reproduced to eight decimal places, namely, we obtain 
—2.2244674718 MeV. As the fitting segment for these calculations, we (arbitrarily) chose 
the interval from Ei = 1 MeV to i?jv = 10 MeV on the real axis. The fitting points. 




(16) 




En — El 




(17) 
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were evenly distributed over this interval and the corresponding momenta were taken as 

kn = \j2mEn/'h^ ■ Table [U shows the meaningful as well as all spurious poles up to = 5. 
It is amazing, but even with only one fitting point we already get a reasonable value for 
the binding energy. 

The second testing potential also describes the 5'-wave proton-neutron interaction, but 
in the singlet state (total spin = 0), 

vSl{r) = -W^'^e^^{-r/R,) (18) 

with = 104.20 MeV and Rq = 0.73 fm. It has a weaker attraction and thus, instead 
of a bound state, it generates a virtual state at E = —0.0660644719 MeV. Similarly 
to the triplet case, this is the only spectral point in the physically reasonable domain. 
The convergence here is even more faster (see Table [21 where we show the results up to 
= 5). With the same choice of the fitting interval and fitting points as we used for the 
first potential, ten decimal places of the virtual state energy are reproduced already at 
N = 7. 

The next example is an exponential hump, shown in Fig. |2j It is positive everywhere 
and therefore can only support resonant states. In the units such that h = m = 1, the 
functional form of this potential is 

V{r) = 7.5r2exp(-r) . (19) 

The exact S-matrix corresponding to this potential, has an infinite number of poles form- 
ing a string that goes down the fc-plane to infinity. The exact locations of the first nine 
of these poles are given in Table |3l These values of the resonance energies were obtained 
using a very accurate method, which is based on a combination of the complex rotation 
and a direct calculation of the Jost function (see Ref. f4j). 

It is naturally to expect that the approximate S'-matrix will reproduce the beginning 
of the string of resonance poles, i.e. the most significant poles that are close to the real 
axis and thus to the fitting segment. This is indeed the case as is seen in Figs. |3] and 
m The first of these two figures shows the distribution of the exact and approximate 
S'-matrix poles over the fc-plane, with the number of fitting points N = 20. The second 
figure shows the same, but for = 30 (few spurious poles that are too far away from 
the origin are not shown). In both cases the fitting points were uniformly placed on the 
interval 1 < < 10, as is given by Eq. (fTTj) . In Figs. [3] and HI the corresponding 
points on the real fc-axis are indicated by vertical bars. As is seen, the distances between 
neighboring bars are decreasing towards the right end of the interval. This is because the 
relation between k and E is not linear. 



9 



The poles of the exact S'-matrix are symmetric on the k-plane relative to the imaginary 
axis (see Fig. [T]). This is a consequence of the Schwartz reflection principle and of the 
fact that the incoming and outgoing spherical waves swap their roles under the operation 
of complex conjugation (for proof see, for example, Ref. |12])- The distribution of the 
Fade poles, shown in Figs. [3] and HI is almost symmetrical with respect to the imaginary 
/c-axis, despite the fact that we fitted the S'-matrix only on a short interval on the right 
hand side of the real axis. This clearly shows that the approximate S'-matrix when 
fitted to the exact points on that interval, possesses the correct symmetry properties. 

Comparing Figs. [3] and HJ we see that, when is increased, the meaningful poles 
converge to the exact ones while the spurious poles change their positions. It is not 
difficult to design an algorithm for automatically selecting the meaningful poles. Indeed, 
as is seen, the meaningful are those poles that are located in the fourth quadrant of the 
fc-plane and do not have a symmetric (or nearly-symmetric) partner at —A; (see Eq. f|T5|) 
and the associated reasoning). 

It should be noted that we did not make any effort to choose the best fitting interval 
and the best distribution of the fitting points. This is a model case, therefore, when doing 
the calculations, we already knew the exact positions of the poles and could distribute the 
points in such a way that the convergence would be much faster. In realistic problems, 
however, no prior knowledge about resonance energies is available. This is why we chose 
the fitting points in a sense "arbitrarily". Despite this, the first few resonances were 
reproduced to a high accuracy. In Table HJ the exact and approximate S'-matrix poles are 
compared for the first five resonances generated by the potential 0191) . 

When applying this method to an unknown potential, one should do it in a few it- 
erations. At the first step the fitting interval and points are chosen arbitrarily and the 
calculations are repeated with different N. Looking at the results, it is easy to identify 
possible meaningful poles. If an exact pole is very close to the real axis (extremely narrow 
resonance) then the corresponding approximate pole at this stage may be found slightly 
above the real axis. So, the poles that are very close to the real axis but are located above 
it, should not be excluded at the first iteration. 

At the second step, one should choose a fitting interval that covers all "suspected" 
resonances. If there is a possible narrow resonance, the density of the fitting points should 
be higher next to it. Repeating the calculations with the adjusted fitting interval, one 
should be able to obtain the resonance energies and widths of at least few most significant 
resonances. 

A good guidance for choosing the fitting interval and fitting points is the energy (or 
momentum) dependence of the scattering phase shift 6e{k), if available. Indeed, we are 
trying to construct such Se{k) that would be as close as possible to the exact function 
Si{k) = exp[2i5i{k)] on the real axis. So, the fitting interval should cover all significant 



10 



variations of Se{k). The density of the fitting points should be high in places where the 
function 6£{k) changes very rapidly, i.e. near narrow resonances. 

The S-wave phase shift function 6o{k) for the potential f|T9|) together with our "ar- 
bitrary" choice of the fitting points ki, k2, ■ ■ ■ , are shown in Fig. [51 As is seen, our 
choice of these points is not the best. They should have been more dense where 6o{k) 
jumps in vr due to the first resonance. We however deliberately stick to the simple set 
of points evenly distributed over the energy interval 1 < E < 10. This is done for the 
purpose of testing the robustness of the method with a poor choice of the fitting points. 
The results obtained show that even with such a choice the most significant resonances 
are reproduced rather well. The quality of the approximation of the exact ^'-matrix on 
the fitting interval with = 30 can be deduced from Fig. O 

It may seem that one can avoid all this hassle with choosing the fitting points by 
simply increasing A^. Such an approach, however, is not viable. The problem is that with 
large A^ the linear system ffTTl) becomes ill-conditioned because the neighboring points kn 
and kn+i are too close to each other and the corresponding equations of the linear system 
are only slightly different. This is similar to a general and unavoidable problem of fitting 
functions by high-order polynomials, where one has to solve the so called Vandermonde 
system, which is ill-conditioned due to the same reason (see Chapter 3 of Ref. [5]). For 
the potential (fT9l) . we found that the instabihties caused by this problem start to show 
up when there are more than 70 points placed on the interval 1 < < 10. 

6 Conclusion 

The proposed method is based on a rigorous mathematical fact that allows us to analyti- 
cally continue the S'-matrix, known on the real axis, to the domains of complex momentum 
where it may have poles corresponding to the bound and resonant states, which can thus 
be easily located. Since the only input information we need, is a table of the S'-matrix 
values along the real fc-axis, the method is universal, i.e. independent of the nature of the 
underlying interaction and the way the table is obtained. 

The potential can be non-analytic that does not allow using the complex rotation 
methods. It can be non-local or energy-dependent, which also makes it extremely difficult 
to apply the rotation. And in principle the S-matrix table can be obtained from experi- 
mental data by means of the phase-shift analysis of the cross section. In all these cases 
the spectral points can be immediately located as soon as the S'-matrix is calculated for 
real k. 

The numerical examples show that the proposed method is stable and accurate. With 
just few matching points, it reproduces the bound states and the most significant reso- 
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nances to the accuracy that is sufficient for any practical purposes. 
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Figure 1: Symmetry properties of the Jost functions on the fc-plane. The dashed lines connect 
the points at which the values indicated next to them are identical. 
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Figure 2: Testing potential (IT9i) in the arbitrary units such that h = m = 1. 
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Figure 3: The exact positions of the S-wave resonance poles (dots) on the momentum plane 
for the potential (HM . and the corresponding poles of the Pade approximation (open circles) 
obtained with 20 fitting points evenly distributed over the energy interval 1 < < 10 on the 
real E-ax\s. The corresponding fitting points on the A;-axis are indicated by vertical bars. 
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Figure 4: The exact positions of the S-wave resonance poles (dots) on the momentum plane 
for the potential ( ITOll . and the corresponding poles of the Pade approximation (open circles) 
obtained with 30 fitting points evenly distributed over the energy interval 1 < < 10 on the 
real E-ax\s. The corresponding fitting points on the k-ax\s are indicated by vertical bars. 
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Figure 5: The exact S-wave phase-shift for the potential ( IT9f ) as a function of the collision 
momentum. The vertical bars on the fc-axis show the points kn, n = 1,2, N {N = 30) 
were the exact and the approximate ^-matrices coincide (the Pade fitting points). 
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Figure 6: Error in the Pade approximation of the S'-matrix for the potential (HM with 30 
points at which the error is zero. 
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N 


Rek (im ^) 


mvk (im ) 


ReE (MeV) 


imh (MeV) 


1 


0.2 X 10-16 


0.1595408980 


-1.0555482734 


0.3 X 10-15 


2 


-0.6 X 10^16 


0.9148722594 


-34.7100271832 


-0.5 X 10-1^ 


0.4 X 10-16 


0.2251639051 


-2.1024785800 


0.7 X 10-15 


3 


-0.1 X 10-12 


-6.5048412082 


-1754.71841570 


0.7 X 10-1° 


-0.9 X 10-16 


0.2305398780 


-2.2040739084 


-0.2 X 10-1^ 


0.3 X 10-1^ 


0.7544907890 


-23.6070608641 


0.2 X 10-12 


4 


-0.9 X 10-12 


1.4615496806 


-88.5852061340 


-0.0000000001 


-0.2 X 10-12 


0.6623217412 


-18.1916485831 


-0.9 X 10-11 


0.2 X 10-1^ 


0.2315901300 


-2.2242014946 


0.4 X 10-13 


0.8 X 10-12 


-0.9942646075 


-40.9956706906 


-0.6 X 10-1° 


5 


-0.2 X 10-1° 


1.3860254952 


-79.6666349435 


-0.0000000021 


-0.1 X 10-11 


0.6586014432 


-17.9878555574 


-0.6 X 10-1° 


0.4 X 10-1^ 


0.2316032946 


-2.2244543701 


0.8 X 10-13 


0.8 X 10-11 


-0.9618425253 


-38.3655990692 


-0.0000000007 


0.0000000207 


-85.5633495380 


-303605.468939 


-0.0001469872 



Table 1: Poles of the approximate 5- wave function 5*0 (A;), found for the potential f|T6|) with 
fitting points ki, k2, . . . , kj^ evenly distributed over the interval 1 MeV < -E < 10 MeV. 
For = 1, the single fitting point corresponds to E = 1 MeV. The boldfaced numbers 
converge to the exact bound-state energy —2.2244674752 MeV. 
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Refc (im ^) 


Im k (im ) 


Rei? (MeV) 


\m.hj (MeV) 


1 


-0.2 X 10-^*^ 


-0.0753661730 


-0.2355520897 


0.1 X 10-15 


2 


-0.8 X 10"^^ 


0.7896265828 


-25.8569655146 


-0.5 X 10-1^ 


-0.1 X 10"^^ 


-0.0409220958 


-0.0694464054 


0.4 X 10-1*^ 


3 


-0.2 X 10-^^ 


0.7187567981 


-21.4238720550 


-0.1 X 10-12 


0.2 X 10-16 


-0.0399961842 


-0.0663393412 


-0.8 X 10-16 


0.4 X 10-12 


-11.9310494904 


-5903.25209234 


-0.0000000004 


4 


-0.9 X 10-12 


-1.1951630087 


-59.2363541815 


0.8 X 10-1° 


-0.2 X 10-15 


-0.0399135508 


-0.0660655059 


0.7 X 10-15 


0.6 X 10-13 


0.6857755943 


-19.5028502331 


0.4 X 10-11 


0.1 X 10-11 


1.4486760697 


-87.0315278526 


0.0000000001 


5 


-0.0000005377 


-259.016625071 


-2782206.21227 


0.0115523186 


-0.5 X 10-1° 


-1.1697429680 


-56.7433434047 


0.0000000053 


-0.7 X 10-15 


-0.0399132475 


-0.0660645021 


0.2 X 10-1^ 


0.2 X 10-11 


0.6850556465 


-19.4619223746 


0.9 X 10-1° 


0.9 X 10-1° 


1.4068361606 


-82.0769256396 


0.0000000102 



Table 2: Poles of the approximate S-wave function Si^ik)^ found for the potential f|T8|) with 
fitting points ki, k2, . . . , kj^ evenly distributed over the interval 1 MeV < -E < 10 MeV. 
For = 1, the single fitting point corresponds to E = 1 MeV. The boldfaced numbers 
converge to the exact virtual-state energy —0.0660644719 MeV. 
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Rek 


Ini k 




r 


2.6177861703 


-0.0048798793 


3.4263903101 


0.0255489612 


3.1300424437 


-0.3571442525 


4.8348068411 


2.2357533377 


3.3983924252 


-0.9972518977 


5.2772798640 


6.7781065905 


3.5914630921 


-1.6639555063 


5.0649296074 


11.9520695757 


3.7383048831 


-2.3317810362 


4.2688602993 


17.4338168679 


3.8534573944 


-2.9922517758 


2.9477816003 


23.0610294625 


3.9448589582 


-3.6424641005 


1.1471837383 


28.7380142741 


4.0173695706 


-4.2816627873 


-1.0966889789 


34.4020435866 


4.0742577101 


-4.9099759809 


-3.7541441225 


40.0090149934 



Table 3: The S'-wave resonance points, A;^/(2m) = — iT/2, of the exact S-matrix for the 
potential iHM . All the values are given in the arbitrary units such that h = m = 1. They were 
obtained using the Jost function method described in Ref. [4^. 





Re A; 


Im k 


exact 

approximate 


2.6177861703 
2.6177861702 


-0.0048798793 
-0.0048798793 


exact 

approximate 


3.1300424437 
3.1300424420 


-0.3571442525 
-0.3571442539 


exact 

approximate 


3.3983924252 
3.3984191164 


-0.9972518977 
-0.9972641751 


exact 

approximate 


3.5914630921 
3.6001916924 


-1.6639555063 
-1.6689691843 


exact 

approximate 


3.7383048831 
3.8771630790 


-2.3317810362 
-2.5129915936 



Table 4: Comparison of the first five resonance points for the potential (IT9l) . obtained using 
the rigorous Jost function method [4J (exact) and the Pade approximation with the number 
of fitting points = 30 evenly distributed over the interval 1 < E < 10 (the units are such 
that h = m = 1). 



20 



